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The electron-electron pair distribution functions (PDF) of the 2-D electron fluid (2DEF) in the 
quantum regime (at T=0) are calculated using a classical-map-hyper-netted-chain (CHNC) scheme 
and compared with currently available Quantum Monte-Carlo (QMC) simulations in the coupling 
range r 3 — l to 50. We iteratively extract the bridge function of the "equivalent" classical 2-D liquid 
in the quantum regime. These bridge functions B(r) are relatively insensitive to spin-polarization 
effects. The structure of the bridge functions changes significantly for r 3 > 6, suggesting the onset 
of strongly correlated clusters. The new B(r), appropriate for the long-range Coulomb potential, 
can be used to replace the hard-sphere B(r) previously used in these calculations. They provide 
accurate classical representations of the QMC-PDFs even at very strong coupling, and probably at 
finite-T near T = 0. 
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I. INTRODUCTION 

The pair-distribution functions (PDFs) of strongly- 
coupled electron fluids contain all the physical informa- 
tion associated with the ground-state static properties 
of such systems. Exchange-correlation energies, phase- 
transitions, and Fermi-liquid parameters like the effec- 
tive mass m* , and the spin-susceptibility enhancement 
(g*) can all be evaluated from the PDFs, as discussed 
below. The static local-field corrections to the response 
functions can also be addressed via these PDFs. 

The PDFs are usually determined by quantum Monte- 
Carlo (QMC) simulations, since standard many-body 
methods become unreliable for densities where the 
electron-sphere radius r s exceeds unity. The r s param- 
eter is also the ratio of the mean Coulomb energy and 
the Fermi energy, and hence is a measure of the coupling 
strength. QMC simulations for the 2D electron system 
were first published by Tanatar and Ceperley-^, and more 
recently by Attaccalite et ali^, and by Drummond and 
Needs^. A transition from the paramagnetic state to the 
ferromagnetic phase was predicted to occur at r s ~26 by 
Attaccalite et al., while weak-coupling theories (Hartree- 
Fock, RPA) predicted such transitions already at very 
low values of r s . In contrast, Tanatar and Ceperley, as 
well as the most recent QMC work by Drummond et al., 
find no such phase transition, where the ferromagnetic 
state is very close in energy, but the paramagnetic state 
remains the ground state. Indeed, direct comparisons of 
the QMC-PDFs of Attaccalite et al., with the more ac- 
curate PDFs of Drummond et al., show slight differences 
which show the need to be cautious about undue claims 
of final accuracy. 

QMC methods have also been used to calculate Fermi- 
liquid parameters^ like m* and g*, but these, and es- 
pecially the m* calculation^ may require further effort 
before a consensus is reached. 

Besides QMC, one other method 5 available for the 
calculation of PDFs of quantum (e.g., electron) fluids 



at arbitrary coupling, temperature and spin polariza- 
tion is based on considering a classical charged fluid at 
an assigned classical- fluid temperature T c f, selected to 
reproduce the correct correlation energy of the quan- 
tum fluid which has two spin species a, f3. The non- 
interacting PDFs 5° g (r) of the classical charged fluid 
are formulated to agree with the analytically known non- 
interacting g° g (r) of the quantum fluid. This is done 
by introducing an effective potential known as the Pauli 
exclusion potential P a ,/3 to exactly reproduce the Fermi 
hole in the parallel-spin g°(r). Once an "equivalent" 
classical map of the quantum fluid is constructed, its 
PDFs are obtained using a classical integral equation, 
namely, the hyper-netted-chain (HNC) equation modified 
to include bridge corrections. This classical-map hyper- 
netted chain procedure was called the CHNC, and we 
showed that the method was surprisingly accurate. The 
classical-fluid temperature of the quantum fluid at T = 
was called the "quantum temperature" T q . In fact, us- 
ing only the correlation energies tabulated by Tanatar 
and Ceperley as the inputs, we constructed the effective 
quantum temperature T q as a function of r s , and showed 
that accurate PDFs of the quantum fluid could be cal- 
culated at arbitrary r s via the classical map 5 . Balutay 
and Tanatar— also presented a closely similar tempera- 
ture map. The advantage of the CHNC approach is that 
it is numerically very simple, and makes many proper- 
ties (spin-dependent properties, Fermi-liquid parameters, 
local-field corrections, finite-T results, etc.) easily com- 
putable with negligible effort. Thus the product m*g* 
can be evaluated^ from the second derivative of the cor- 
relation energy with respect to the the spin-polarization 
£. The effective mass m* can be obtained^ from the sec- 
ond density-derivative of the free energy which can be 
evaluated from the temperature-dependent PDFs. 

However, the classical map of the 2D- electron fluid 
leads to a Coulomb fluid for which the simple HNC is in- 
adequate. The irreducible three-body and higher-order 
contributions which are lumped into the bridge term 



B at p(r) are very difficult to calculate directly. The need 
for bridge functions appears in many areas in the theory 
of Coulomb systems^. Hence, as is customary^ in our 
earlier work we used a hard-disk model, based on solu- 
tions of the Percus-Yevik equation 10 . However, the avail- 
ability of an extensive set of QMC-generated PDFs for 
the 2D-electron fluid presents the possibility of extracting 
accurate Coulomb- adapted bridge functions B(r). These 
are implicit functions of the PDFs themselves^. Thus, 
in this paper we present an iterative procedure for ex- 
tracting the 2D-bridge functions appropriate to a clas- 
sical Coulomb fluid at a given r s , and at the effective 
classical-fluid temperature T c y, and yielding the quan- 
tum PDF at the given r s and at T = 0. 

II. HNC AND CHNC METHODS 

The HNC equation^ and its generalizations, coupled 
with the Ornstein-Zernike equation have lead to very ac- 
curate results for classical charged-particle interactions. 
In the following we use indices i,j which could be spin 
indices or other species-identifying indices. The exact 
equations for the PDFs are of the form: 

gij (r) = e-A/^W+iVyM+ByM (1) 

Nij = h l3 {r) - Cij (r) (2) 

M r ) = 9ij(r) - 1 (3) 

Here </>ij(r) is the pair potential between the species 
and Nij (r) is the nodal function, while cy (r) is the direct 
correlation function connected to hij by the Ornstein- 
Zernike (O-Z) equation. The pair potential <fiij( r ) is the 
the sum of a diffraction corrected Coulomb potential and 
the Pauli exclusion potential (see Eq. 1-2 , ie&). 

If the bridge function Bij(r) were set to zero we have 
the HNC approximation, adequate for systems where 
the kinetic energy dominates strongly over the poten- 
tial energy. The bridge function brings in many-body 
cluster interactions beyond the diagrams of the hyper- 
netted chain expansion. A systematic investigation of 
the 3-D one-component plasma was given many years 
ago by Lado et al 9 . As the bridge interactions involve 
many-internal interactions (within the cluster) averaged 
over, they are somewhat insensitive to the exact form 
of the pair-interactions, and their spin-states. Thus it 
was showni^ that the analytically available Percus-Yevik 
hard-sphere bridge function could be used for most flu- 
ids, to closely reproduce the Monte-Carlo PDFs available 
at the time. The hard-sphere radius of the model bridge 
function was in effect an optimization parameter in such 
approaches to liquid structure. 

The PDFs of the classical 2D electron fluid became 
relevant to 2D quantum fluids when Laughlin^ intro- 
duces his plasma mapping of 2D fractional-quantum Hall 
(FQH) fluids to classical plasmas. The use of a suitable 
bridge function was found to be essential if high accu- 
racy was to be obtained via the plasma mapi^ for the 




FIG. 1: (Online color) The QMC pair-distribution function of 
a fully spin polarized = 1) electron fluid at r s = 10, T = 
is compared with those calculated from CHNC using a hard- 
sphere bridge function and with no bridge function what so 
over 



hierarchy of FQH states. The temperature of the classi- 
cal fluid in Laughlin's classical map is directly related to 
the filling factor, via the form of the many-body wave- 
function proposed by Laughlin (Laughlin used a simple 
HNC without bridge corrections). There is no external 
magnetic field in the 2D systems studied here. It turns 
out that the temperature of the classical charged fluid 
whose PDF agrees with the electron (quantum) fluid at 
T = can be expressed^ as a function of the density 
parameter r s . 

We founds that the use of hard-sphere model bridge 
functions gave good agreement around the first peak of 
the PDFs, while the more distant oscillations were weaker 
than those in the QMC-PDFs. The PDF of the fully 
spin-polarized T = electron fluid at r s = 10, calculated 
using CHNC without bridge, and CHNC with a hard- 
sphere bridge function^, and the bench-mark QMC-PDF 
are shown in Fig. [TJ 

The bridge function B(r) is a function of the PDFs 
themselves, and hence their extraction from QMC simu- 
lation results is by no means obvious. However, given a 
choice for the effective pair potential and classical-fluid 
temperature implicit in the CHNC /30(r) y , and the as- 
sumed applicability of the modificd-HNC equation and 
the O-Z equation, a bridge function can be evaluated 
from the QMC-j(r). This evaluation requires separating 
out the long-range tails of the pair-potentials, direct cor- 
relations etc., in r- and k- space, using numerical Fourier 
transforms for the short range parts, analytical formu- 
lae for the long ranges parts and reassembling the results 
in r- and fc-space. Since these mathematical machinery 
are already available in the algorithms for the CHNC, we 
present here a simple iterative scheme equivalent to the 
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above, but using the CHNC code itself. 

The B(r) functions obtained from such a procedure, 
together with the /3<j>(r) of CHNC present a complete, 
accurate, classical representation of the quantum PDFs 
obtained by QMC. It is hoped that such bridge functions 
and classical maps can then be used in regimes where 
QMC simulations are not possible or easily available, as 
discussed below. 



III. EXTRACTION OF THE BRIDGE 
FUNCTION FROM QMC DATA 

Given a density n = l/far^) specified by an r s value, 
the target QMC-PDFs that we use are those of Attac- 
calite et al., as parametrized by Gori-Giorgi et alJ£. We 
drop the species indices i,j except when required, and 
indicate the PDFs from QMC and CHNC by g(r) and 
gchndr). Given the hypothesis that the QMC-PDFs can 
be represented by classical forms, a nodal function and 
bridge function corresponding to the given target g(r) 
should exist. Thus the QMC and the CHNC PDFs sat- 
isfy 

log [g(r)] = -P<Kr) + N(r) + B(r) (4) 

Both N{r) and B{r) are implicit functions of g(r). How- 
ever, all the terms in the CHNC form 

log [g c hnc{r)\ = -P<p(r) + N chnc (r) + B chnc {r) (5) 

are known. Also, it has been found from previous com- 
parisons^ of g(r) and g c hnc(r) that they agree closely, 
even when B c ? mc (r) was taken from a hard-disk model. 
Hence we assume that 2V(r) can be replaced by N c h n c as 
a first approximation. Then we easily obtain an initial 
estimate of the bridge function contained in the QMC- 
PDF. 

B{r) = log [g(r)} - log [g c h n c(r)} + B chnc (r) (6) 

Thus, starting from the hard-sphere model of B c h nc (r) we 
obtain a a standard mixing procedure to construct a new 
Bchncif), and hence a new g c hnc, and so on. In the small- 
r region the value of g(r) becomes negligibly small and 
hence the extraction of the difference between two log- 
arithms becomes numerically unsatisfactory. However, 
as may be surmised from Eq. [Bl the calculated g(r) are 
found to be insensitive to the form of B(r) for r < r s . 
Thus even a simple polynomial extrapolation (connecting 
the r/r s < 1 region with, say the region 1 < r/r s < 1.5) 
or the hard-sphere model itself may be used. The itera- 
tive procedure is even insensitive to slight discontinuities 
at the connection point (although of course, discontinu- 
ities should be avoided). These bridge functions will be 
called "Coulomb bridge functions", to distinguish them 
from the hard-sphere bridge functions BHs( r )- 

The extraction procedure for the Coulomb bridge func- 
tions Bij(r) from the QMC-PDFs, i.e., gij(r) is found to 



be very efficient, and the iterative inclusion of the ex- 
tracted bridge functions leads to rapid convergence in 
reproducing the target gij(r). It is easiest to consider a 
fully spin-polarized (i.e., £ = 1) electron fluid as it is a 
one- component system, with just one bridge function. 
A set of bridge functions for a range of r s values which 
reproduce the QMC pair-distribution functions at £ = 1 
when used in CHNC are shown in Fig.[5J It is found that 
as r s increases from 0.3 (not shown) to unity, the B(r) 
remains more or less unchanged. From then onwards, 
esp. after r s — 3, the development of the first peak in 
the PDF produces an oscillatory structure in B(r). This 
trend continues till about r s — 6. However, now the 
deepening of the second trough in the PDF leads to a 
complete qualitative change in the the form of the bridge 
function, as seen from the panel (b) of Fig. [5J The new 
deep trough seen near x = r/r s ~ 3 is already there as a 
weak trough in the moderately coupled fluids of panel (a) . 
In panel (a), the deep trough is in the r/r s < 1 region and 
it is not of much importance in the HNC equation. In 
panel (b), when the deep trough near x ~ 3 develops, the 
small x region rises rather steeply. But this rise has only 
a small numerical significance in the HNC, since the pair 
potential and the Pauli potential become very large and 
dominant for small x. In contrast, the structure in B(r) 
for larger x becomes significant, as the pair-potential ef- 
fects fall off with increasing x. These results show that 
the short-range structure of the fluid undergoes signifi- 
cant and subtle changes in the regime r s > 6. Previous 
studies of the local-field corrections to the response func- 
tion of the 2D electron liquid had unraveled interesting 
characteristic o 16 ' 17 which may well be related to the onset 
of strongly correlated clusters for r s > 6. 

Another system which effectively reduces to a one- 
component system is the paramagnetic electron fluid 
(C = 0) with equal amounts of up-spin, and down-spin 
species. As a test of spin insensitivity, we can use the 
bridge functions determined from the £ = 1 case for the 
C = case and see if the CHNC-PDFs reproduce the 
QMC-PDFs. This is in fact very nearly the case, show- 
ing that the spin-polarization dependence of B(r) is very 
small. This also shows that, if desired, one may intro- 
duce just one bridge function £?(r, £) for all i,j com- 
ponents, where the latter is based on linearly interpo- 
lating between B(r,£ — 1) and B(r,£ = 0), since they 
are very similar. The differences are found in the near 
r ~ r s region close to the first peak, where the Pauli- 
exclusion effects, diffraction effects etc., are comparable 
to the Coulomb repulsion effects (see the discussion of 
Fig.®. 



IV. DISCUSSION 

Currently, the most sophisticated QMC calculations 
for the quantum 2D electron liquid are those of Drum- 
mond et al. In practice, the differences between those 
of Drummond et al., and those of Attaccalite et al., are 
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FIG. 2: (Online color)The 2D Coulomb bridge functions 
needed for a classical description of the pair-distribution func- 
tion of a fully spin-polarized (( = 1) 2D electron fluid, for 
r a = 1 to 50. Panel (a) shows the weak to moderately coupled 
regime. Panel (b) shows that the nature of the B(r) changes 
rapidly near r a — 6 and r s — 7, when the second trough in 
the PDF begins to develop. The PDFs are insensitive to the 
behaviour of these functions for r/r a <~ 1.3. 



FIG. 3: (Online color) A comparison of the 2D electron-liquid 
PDFs of Drummond et al^ with those obtained from the 
CHNC, using the bridge functions iterated from Eq. In 
the case of r 3 — 30, the Drummond et al. data are for ( — 0, 
and the effect of using the £ = 1 bridge function on the PDF 
(dashes) in CHNC instead of the correct B(r) at £ = is 
also shown (PDF, Boxes). The PDFs at r 3 = 30 have been 
displaced upwards for clarity. 



irrelevant except for dealing with very small energy dif- 
ferences indicative of a phase transition from the para- 
magnetic to the ferromagnetic state. The presence or 
absence of a ferromagnetic transition would also modify 
the predictions regarding the enhancement of g*. The 
lack of such a transition in the results of Drummond et 
al., further illustrate the difficulties inherent in assuming 
that a given set of QMC data are essentially the defini- 
tive result. Given the very slight energy differences be- 
tween the two phases, any technical improvements can 
modify the conclusions. Since CHNC is based on us- 
ing some of the QMC data as inputs, and then applying 
CHNC to other instances, the conclusions of CHNC also 
get modified as the input QMC data are modified. Thus, 
in our study of the two- valley 2D electron gasi& found 
in Si metal-oxide field-effect transistors (MOSFETS), we 
used hard-sphere bridge functions, and the classical-fluid 
temperatures based on the QMC correlation energies of 
Tanatar and Ceperley as inputs to the CHNC. It was 
found that the two- valley systems becomes critical at suf- 
ficiently high r s values, when the effective mass and the 
spin susceptibility become very large. This property is no 
longer found to be the case when more accurate bridge 
functions are used. 

The pair distribution functions calculated from CHNC 
using the bridge functions extracted from the Gori-Giorgi 
et al. fit to the Attaccalite g(r) are, of course, essen- 
tially identical to the target g{r). Thus a comparison of 
CHNC-g(r) using the fitted B(r) is a check on our fit- 
ting process as well as an indirect comparison of Refs. 2 
anc£. It also enables us to check on the applicability 
of the B{r) obtained from the polarized-2D system for 



predictions on the unpolarized system. In fig. [3] we com- 
pare CHNC- g(r) with numerical data from Drummond 
et alJ£, for r s = 20, C = 1, and r s = 30, ( = 0. In the lat- 
ter case, we have also calculated the CHNC PDFs using 
the appropriate bridge function B(r, ( = 0), and also the 
inappropriate B(r, C = 1) to display the relative insensi- 
tivity of the PDF to any spin-polarization dependence in 
B(r). 

It should also be noted that we have not modified the 
quantum temperature T q of the classical map given in 
Ref£ in extracting the bridge functions. In the classi- 
cal map of Balutay et al., a pure-HNC procedure and 
a different T q were used. If it were extended to include 
the bridge terms, qualitatively very similar results are 
obtained for the bridge functions which are now slightly 
different functions. 

Once a family of bridge functions for each r s and T q is 
obtained, the g(r) for each r s as well as all the smaller 
r s data are used in the adiabatic-connection formula for 
the correlation energy E c . Then the QMC value of E c 
will be accurately recovered, and hence we need not re- 
vise the preset quantum temperature T q that was used 
in the classical map. That is, given a T q of a classical 
map of the CHNC type, the use of the Coulomb bridge 
function based on that T q instead of the original bridge 
functions used when forming the T q , does not require any 
refitting the T q . This is a very convenient conservation 
property for implementing CHNC calculations with im- 
proved bridge functions. 

The advantage of the CHNC procedure over those 
of QMC is its simplicity of implementation, as well 
as its easy applicability to finite-T, finite-^ situations. 
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Thus CHNC studies of multi- valley systems (e.g., Si- 
MOSFETs), nanostructures etc., can be attempted. 
Finite- T studies become possible if the B(r) remain valid 
even for fmite-T, at least in the near T — region. An es- 
timate of the transferability of the T — bridge functions 
to fmite-T may be obtained by noting that temperature 
acts to level out the oscillations in the PDFs, while the 
B(r) attempts to sharpen the oscillations. 

Some years ago we showed^ that very interesting spin- 
dependent effects arise in the T = to T — E F region 
due to the interplay of spin, onset of partial degeneracy, 
Coulomb correlations and temperature effects. Those 
studies were carried out with the hard-sphere bridge 
functions available at that time. Such temperature- 
dependent studies are still beyond the reach of QMG 
methods. The accuracy of those earlier f?(r)#s-based 



studies at fmite-T can now be further examined using 
the new bridge functions. 

Finite temperature studies of the free energy of the 2D 
function of density and spin polarization can also 
be used to determine Fermi-liquid parameters like the ef- 
fective mas o 4 i 21 and the spin susceptibility. In CHNC we 
directly determine m* from the second density derivative 
of the excess free energy with respect to the density^. 
Thus, the availability of these new Coulomb bridge func- 
tions which accurately describe the T = PDFs provide 
greater confidence in the results that may be obtained 
for the Fermi-liquid parameters in future studies. 

In conclusion, we have presented an accurate classical 
representation of the pair-distribution functions of the 
quantum 2D electron fluid at T = 0, for arbitrarily strong 
coupling and arbitrary spin polarizations. 
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